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Recent analysis of spontaneous breaking of SU(2) as a continuous gauged flavor symmetry 
to non-abelian discrete group as a residual symmetry motivates its consideration as 
a group for stabilizing dark matter. We determine the region of hidden sector particles 
compatible with relic dark matter density and present prospects in which these particles 
could be observed using direct and indirect dark matter search experiments. 

I. INTRODUCTION 

Closeness of Gauge-couplings at the Unification scale, neutrino masses, structure of lepton mix- 
ing matrix and hierarchy of quark mixing angles, dark matter, dark energy and baryon asymmetry 
of the Universe are some of the still unexplained experimental facts. While framework of solving 
them all doesn't exist yet, we can pick one issue and offer possible solution which may lead towards 
solving more of them. This can be done on different levels: choosing new fields and adding them 
within appropriate representation, extending the gauge symmetry i.e. following the gauge route, 
adding supersymmetry or extending and modifying basic concepts like quantum fields and nature 
of space-time, increasing levels of speculation with each step we take. Here we will focus on the 
search for the dark matter (DM). 

There are many evidences of DM existence \T] and reasons for expecting it to be of particle 
nature. Such a particle would have to be neutral, stable over cosmological scales and provide 
current DM relic density. It is reasonable to assume that the group used for accommodation of 
DM and SM extension simultaneously describes more observations. Special structures of the lepton 
mixing matrix and the quark mixing matrix give quite a strong hint for a flavor symmetry G broken 
in a non-trivial way. Requirement that at least two fermion generations can be unified by G leads 
us to the conclusion that the best choice for G is a discrete, non-abelian symmetry. [2j 
Since the discrete symmetry that stabilizes dark matter must be a remnant of a spontaneously 
broken gauge symmetry [4J, in most cases the used symmetry is U{l)x gauge symmetry which is 
broken to discrete Z2 subgroup. Preservation of Z2 symmetry stabilizes hidden sector in which 
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the lightest particle is usually DM candidate. Examples of Z2 symmetry are R-parity in Minimal 
Supersymmetric Standard Model and T-parity in Little Higgs Model [5]. 

Simplest way to analyze extension of SM with particles in hidden sector stabilized with Z2 
symmetry is to add one singlet scalar to SM particles which has been discussed in [6J. 

However, DM can as well consist of more particle types. In this case \J(\)x gauge symmetry 
is broken to group with N > 2 and fields associated to its irreducible representations create 
hidden sector of the model. Stabilization of DM with Z^ symmetry has been studied in [7J where 
the DM candidate was neutral complex scalar singlet, while DM phenomenology with discrete Z2, 
Z3 and Z4 stabilizations of one extra sclalar doublet and one extra complex scalar singlet was 
investigated in [8]. In [9], a non-Abelian discrete group 1)3 was recently employed to stabilize 
hidden sector consisting of one or two component DM. 

Here we investigate the stabilization of one component DM with quaternion group, Q^. It 
has been shown recently that group is the only non-Abelian discrete group which arises from 
spontaneous breaking of continuous (possibly gauged) flavor symmetry SU(2) [13j. It is a double 
valued counterpart of the dihedral group D2, also called and it has been used in series of articles 
by M. Frigerio and collaborators |10| . |12j for new understanding of family symmetries, Vckm 
and Umnsp matrices and correlation of Amf2 to ^23 ■ 

We assign new particles to non-trivial irreducible representations of the group and stabilize the 
hidden sector preserving the symmetry after electroweak symmetry breaking. We add two complex 
scalars in a doublet and three real scalars in singlet representations and we can compare the model 
with one component case of DM stabilization with Da, Z3 and Z2- With that choice of particles 
the most similar model turns out to be the one based on 1)3 symmetry and in this article we follow 
closely their analysis of one-component DM model. In section II we present the group (54, its 
irreducible representations and decomposition of tensor products. Section III introduces potential 
for DM in model, conditions on the model to achieve potential stabilization and new processes 
which take place in the early Universe. It also shows the potential terms and reactions which 
differentiate the from the 1)3 DM stabilization. Considering annihilation reactions, decays 
and DM conversion we use micrOMEGAs package |16|.|18j.[T7j to present its mass dependence on 
the strength of Higgs portal coupling and investigate prospects of direct and indirect detection of 
lightest DM particle in this model. In section IV we conclude. 
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II. Qa group 

Qi contains eight elements a"^h^ for m=0, 1, 2, 3 and k =0, 1, four one dimensional irreducible 

representations l-i , 1 i-, 1 (= = 1,...,4 respectively) and one pseudo-real two 

dimensional representation 2. The generators of one dimensional representations are a = 1, h = 

-1, h = i for I4 and for two dimensional 

, satisfying = 6^, aba = b. Since 

-sts a similarity matrix U between the 
transformation matrix and its complex conjugate, such that 
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The tensor product of two doublets \ ^ \ decomposes in 
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where si^2 G {+, — }) s" = sis\ and = 8282- Dihedral groups and their double covers are nicely 
described in [2], [HJ. 
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III. STABILIZATION OF DARK MATTER 

We want to find the best way to arrange particles in Q4 irreducible representations to differ- 
entiate from DM stabilization with Z2, and groups. It is clear from the multiplication 
rules that Lagrangian extended with particles from singlet representations preserves additional Z2 
symmetry. To address a nontrivial case we add real particles r/j i=2,3,4 in each nontrivial singlet 

representation and complex X = I 1 in a doublet of Q4, which amounts to five new degrees of 

freedom. Scalar sector of the Lagrangian for DM stabilized with Q4 symmetry is invariant under 
Q4 and gauge symmetries and beside SM Higgs field, it has kinetic terms and potential. It can be 
constructed using multiplication rules of singlets and doublet eq.Q-Q and reads 

(XtX)2 + ai3 [m{X^Xh + H.c] 
+a22 [ri2{XX)2 + H.c] + 023% {XX)^ + a2m (XX)^ + a3^7Jf{X^X) + aii7]f{H^H) 
+a^{H^ H){X^ X) + \i{H^HY + \i+i{X^ XUX^ X)i + A5 [{XX^)2{XX)2 + H.c] 

Where an coupling exists for i=2,3 while the term with 014 vanishes. 

. (A 

The symmetry is preserved after electroweak symmetry breaking (H) = ^ I \ Av) — 

v) 

0, (x) = where = —mf/Xi ~ (246 GeV)^ is squared Higgs vacuum expectation value. For 
simplicity we assume CP conservation, so all the couplings in potential are real. Phenomenological 
constraints on the couplings come from demanding that the potential is bounded from below and 
that the electroweak vacuum is the global minimum of the potential. This gives the conditions 
on the terms with quartic couplings which must be positive, while the terms with an, a2i, asi, a^i 
for i = 2, 3, 4 and a^ must be greater than some minimum (negative) value. The parameter /i is 
determined by the condition > 0. 

A model based on the Q4 symmetry has all the terms in the Lagrangian like the model based 
on the Z3 symmetry besides the term + X*"^- The theory based on an Abelian discrete Z3 
symmetry is most similar to the recent one based on the group D3. The minimal theory based 
on the D3 symmetry has nontrivial interaction predicted by the term ^rj (x^ — x*^) which allows 
the decay rj — )• 3x3% • Where rj transforms as a nontrivial singlet and x as a doublet of irreducible 
representations of group. New terms in the theory based on the Q4 symmetry are those with 
coupling an and a2i which allow the decay of r/j particles. For i = 2,3 the decays are into XX 
and x*X* for i = 4 into x*X- There are also new terms that govern reactions which don't 
contribute to change in number and species of DM particles. 
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After expansion v ^ v + h around the vacuum from potential we can read these physical masses: 

ml = 2Xiv'^ 

= + (^) 

To examine the thermal relic abundance of added particles in Q4 model we use the micrOMEGAs 
package |16j.|17j.|18j. Relic abundance of the particles is governed by Boltzmann equations which 
account for the change of the total number of particles we are interested in. 

From five new added particles the lightest one is taken to be our DM candidate. If the decays 
^2,3 — >• XX) ^2,3 — >• X*X*) ^4 ~^ X*X are kinematically allowed this is x particle. If these decays 
were not allowed there would be five DM candidates and calculation of relic density would have 
to be treated using five coupled Boltzmann equations. This number reduces to four because of 
CP invariance of the potential Q = n^* so one needs to solve only Boltzmann equation for 
and total relic density becomes ^x+x* ~ ^^x giving the relic abundance Qh'^ = fi^+^./i^. Under 
certain circumstances when the second-lightest particle is only slightly heavier than the lightest 
one if certain condition (see (6) below) is guaranteed two Boltzmann equations can be reduced 
to obtain one and it is possible to calculate relic density via standard methods [21J. In this case 
there can be presence of co-annihilations if the model allows it. Reactions which contribute to 
Boltzmann equations in calculation of relic density are 

(a) annihilations into SM particles: rjiirji — Xsm, ViVj ~^ ^SM, X*X ~^ Xsm, XX ~^ Xsm 

(b) semi- annihilations: x^2,3 X*H, x^4 — ^ xH, XX ~^ Hr]2,3, X*X ~^ Hrj/^ 

(c) DM conversion: 

(cl) dark annihilations: rjir]i rjjrjj, r]ir]i x*Xi X*X ViVi, X*X ??2^3, 

XX X*X*, XX ??2,3??4, 
(c2) dark co- annihilations: r^x m,3X*, m,3X ViX*, ?72,3^4 XX, V2,3m X*X*, 

mm x*x 

(d) decay: r/2,3 ^ XX, m,3 X*X*, V4. ^ X*X 

The reactions with x analogous to reactions with x* so we list only those for x Vi- 
Quaternion symmetry allows reactions of annihilation into SM particles (a), semi-annihilation (b) 
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recently studied in |21] . DM conversion (c) among which we recognize "dark" co-annihilations and 
"dark" annihilations, and decay (d). In this work we focus on reactions of type (a) with reactions 
(d) necessary included to make sure we have one DM candidate. DM conversion reactions in 
D3 model were shown to change the number of heavier particles while they didn't influence the 
number of lightest DM particle or change the relic density of lightest DM particles for chosen set of 
parameters. The lightest DM particle would be influenced by these type of reactions if the masses 
of the heavier and the lighter particle were nearly degenerate [9]. 

Annihilations. Reactions of annihilation of DM particles into SM particles are common to 
scalar DM models with Higgs portal interaction. In the case m^^,,,. << mn annihilation cross- 
section depends on the decay width of SM Higgs boson in all kinematically allowed final states 
XsM = ZZ,WW,bb... of SM and of the hidden sector ('dark' annihilations). The decaying Higgs 
boson is a virtual particle with a mass of two times the mass of annihilating particles. The cross 

2 

section of annihilating particles in the region niH » m^,r]i reduces to |6J {crv)^^-^^^^ ^^^2 
depending primarily on 05 coupling and DM mass ruy.. To determine the value of Higgs portal 
coupling to X we use the requirement that x gives the right relic abundance within the WMAP 
observations 0.09 < il/i^ < 0.13 like done in [22j while we assume ^72, 3, 4 have already decayed earlier 
to lighter particles. 

The dependence of 05 on for the range of the from 4 GeV to 1000 GeV, often measured 
in experiments searching for DM, is shown on the Figure [T] in the ruy. — 05 plane for mass of the 
rj2 particle > 2TeV so that reactions (d) are allowed for every ruy.. The band shows values for 
which relic abundance is within the WMAP values. In the region under the band 0, > ilwMAP 
and above the band there have to be other reactions strong enough to source the relic abundance 
to give correct 0,. The strength of the coupling as in the region around Higgs resonance niy. w ^ 
is very weak compared to the cases where << and niy. >> and x particles annihilate 
there very efficiently. In this region the value of Q5 does not have to be very strong to account for 
the viable cosmology. 

To analyze dependence of annihilation channel x*X on masses and we observe dependence 
as — for m^-, « ra^ and m^. , » with i=2, 3, 4. To be sure reactions (d) are 
allowed we keep < 2m^. . This choice of parameters for both regions gives good agreement with 
the graph on Figure ([T]). Reason for this is that reactions of the type (a) strongly dominate over 
(c) . For example, cross-section of reactions XX* ~^ "HiVi is of the order of ~ 10"^'^ pb for center of 
mass momentum 500 GeV while dominant SM channel for this set of parameters W'^W~ is ~ 10~^ 
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FIG. 1: Dependence of as coupling on the masses of DM particle to^. Along the band values of relic density 
are within the WMAP allowed region 0.09 < flh^ < 0.13. The region below is excluded > J7 



DM rp, 

WMAP ■ ^ ^1*= 



resonance at sa allows for value to be of order of 10 



pb for m^, m.^. > mn region (for m^=170 GeV and m.^. ~ 400 GeV). 



Dark matter conversion processes are reactions in which number of dark matter particles remains 
the same. There are two types of processes. The "dark" annihilation reactions (cl) are 2 — )• 2 
processes in which the number of DM particles of a particular species changes by two units. The 
reactions (c2) can be regarded as dark co-annihilations. In these reactions the DM particle of one 
species annihilate with DM particle of another species. 

If we had two hidden sectors each belonging to its particle type, allowed decay of the second light- 
est particle to the lightest one would make the lightest one DM candidate forming one-component 
DM. If this decay doesn't exist or it is not kinematically allowed, we have two-component DM 
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and two hidden sectors, with rehc density of each of them obtained by calculating two coupled 
Boltzmann equations. However, if the mass of the second lightest particle becomes smaller than 
1.5 times the mass of the lightest one [TC], |23j the inclusion of co-annihilations lowers the relic 
density and makes possible the simple use of one Boltzmann equation. This is due to possibility 
of introducing the condition 

^ ^ (6) 

for rii number density of dark matter candidate species: i=X5X*5^2 and n^'^ their number density 
in thermal equilibrium with background particle density n = "^^Lini. Increasing the number of 
particle types increases the number of hidden sectors contributing to total relic density and the 
number of coupled Boltzmann equations we need to solve to find individual relic density of each 
sector. If masses of all particle types from hidden sectors are smaller than 1.5 times the mass 
of the lightest particle the situation becomes simple since the number of Boltzmann equations 
reduces to one. In this case i=x, X*5 '^2, ??3, '?/4- It has been shown for a particular model [23] 
that co-annihilations can lower Qh"^ by factor of 6. The condition (6) cannot be guaranteed when 
semi-annihilations are present [21]. This assumption that DM particles which contribute to co- 
annihilation processes are in relative thermal equilibrium in the early universe is also taken in the 
micrOMEGAs calculation. 

To investigate the influence of almost degenerate masses in quaternion model we choose masses 
of DM particles very similar as it is recommended for typical freeze-out temperature T = 
while the recommended mass of the heavier particles, rji for i=2, 3, 4, is m^. < 1.5my. under 
which co-annihilations become included. In the region of lower masses this effect is shown on the 
Figure [2] Since Qh'^ depends on Higgs portal coupling we show the behavior of depending on 
the mass of the second lightest DM particle approaching to mass of DM candidate for couplings 
0:5 = 0.25,0.27, and 0.29, determined from Figure [T] The couplings which start to show effect are 
those with a/^i with i=2,3,4. Dominant reactions from the "dark" sector are those governed by 
Higgs channel. Analysis shows that the effect of couplings 042, 0:43 and 044 is visible for the entire 
range of masses when coupling 05 is smaller than 0.001. 

Direct detection. One of the ways to search for dark matter is measuring the nuclear recoil 
induced by DM scattering off nuclei. For the scalar interactions, when DM candidate couples to SM 
particles via Higgs boson we have spin-independent (SI) interactions. The best nuclei for detecting 
scalar interactions are heavy ones because scalar interactions add coherently in the nucleus |18j . 
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FIG. 2: Dependence of the mass to^^ (second lightest particle) for = 9 GeV on relic density. Approaching 
the mass of rji to the mass of the lightest particle results in contribution of type (c) reactions. 

There are many experiments searching for DM particles in this manner. DAMA [24J has reported 
an excess in annual modulations effect favoring light scalars as DM in agreement with a small excess 
of events by CDMS, EDELWEISS and CoGeNT experiment who have also reported a possible signal 
|25j.|26j. However, the null hypothesis of CoGeNT provides similar reduced chi square as for dark 
matter particle of 9GeV with asi = 6.7- 10~^^ cm^. Two events XENONIOO observed in the region 
for WIMP search can as well be coming form background fluctuations. The probability for this in 
the expected region for WIMP search is 26.4%. They have provided the most stringent limit for 
> 8 GeV with a minimum of o" = 2.0 • 10~^^ cm^ at = 55 GeV [27J. Figure [s] shows the 
obtained spin-independent cross sections for annihilation of x ^i^d x* particles for the case when 
we have annihilation reactions and decays. 
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FIG. 3: The figure shows spin-independent cross sections for reactions of (a) type, with m,^ = 2100 GeV 
chosen to allow reactions (d) for every m^. The cross-sections (black) are obtained varying Higgs portal 
coupling as as shown on Figure [l] to be consistent with WMAP measurements 0.09 < flh^ < 0.13. Full red 
line shows restriction from CDMS and EDELWEISS experimental results [53], the dotted one the results 
from XENONIOO from 225 live days |27, and the green one is region investigated by CoGeNT. 

The parts of parameter space which are allowed under these conditions are heavy DM particles 
> 165 GeV, very light ones and (small) range of masses around Higgs resonance ~ 60 GeV. 
CMS and EDELWEISS experiments restrict the higher masses as the XENONIOO but allow masses 
around the resonance. Inclusion of additional reactions like semi-annihilations would provide a 
wider range of WMAP-restricted values and would possibly decrease the cross-section allowing 
a larger area of the parameter space. Light DM masses favored by CDMS, EDELWEIS, XENONIOO 
and DAMA annual modulation can on the basis of direct detection measurements be accommodated 
in this model. 
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Indirect detection. Primordial reactions which took place in the early Universe take place 
nowdays in regions assumed to consist of dense DM such as galactic center, "galactic halos or 
extra galactic region" . Annihilation of DM in these regions results in a pair of SM particles which 
hadronize and decay into the stable particles and then propagate through Universe. When these 
particles are photons, neutrinos and anti-matter fluxes they can be detected by indirect detection 
experiments. One of the sources observed by Fermi-LAT are dwarf spheroidal galaxies (dSphs) 
of the Milky Way, DM dominated systems, that makes them good targets for DM searches [28] . 
Gamma ray signal from these measurements determines upper limits on DM annihilation cross 
section. In determination of gamma ray flux (in number of photons per cm?s sr) 

(TV 

cp^{E,<t>) = —f^H {((>), (7) 

model dependence enters in DM mass, annihilation cross section av and annihilation rates in 
different annihilation channels. Here /■yiE) = is energy distribution of particle 7 produced 
in one reaction and H{<j)) is integral of the dark matter density squared over the line of sight and 
dependent on the halo profile which is in our case Navarro- Frenk- White (NFW) profile [17j . 

Figure [4] shows the (av) as a function of for masses from 4 GeV to 1000 GeV. For small 
masses < 5 GeV dominant annihilation channel is into cc quarks which for the masses above 
5 GeV become bb. This annihilation channel prevails until around 75 GeV after which x particles 
dominantly annihilate into W~^W~ bosons. We also show the upper limits measured by Fermi- 
LAT observing dSphs with bb and W~^W~ main annihilation channel with later one for heavier 
masses. For the mass range from 95 GeV to 1000 GeV they give upper bound on (av) value from 
the order of ~ 10~^^ to the order of ~ 2 • 10~^^ [28j. The most stringent constraints come from 
DM annihilation measurement from the region of Galactic Center (GC) by Fermi-LAT [29]. They 
exclude the light masses < 20 GeV of the model while the measurements of diffuse gamma ray 
emission from Galactic Halo (masking out Galactic Plane) exclude the masses lighter than around 
15 GeV [30J. We can notice that investigation of both GC and dSphs allows for heavier masses of 
the model. 

Observations of Coma Cluster of Galaxies produced a little less stringent constraints on {crv). 
They are of the order of ~ 10~^^, ~ 10"^*^ [31] in the lowest region of parameter space. 
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FIG. 4: The figure shows annihilation cross section in dependence on for the WMAP allowed values of 
a5 and values of rriri- — 2100 GeV for i=2,3,4. The red solid line shows upper limit for bb and dotted red line 
for W^W^ dominant annihilation channels measured by Fermi-LAT from observing dSphs and obtained 
using combined analysis |28| . Green lines are upper bounds from measuring the gamma rays from Galactic 
Center (GC) for bb and W^W~ shown with solid and dotted line respectively [29], while the pink line is 
obtained from diffuse gamma ray emission from Galactic Halo |30) . 



IV. CONCLUSION 



We have shown that this very simple model based on non-Abelian finite group with one SM Higgs 
doublet can with certain choice of parameters be very similar to D3 model with one component 
DM. From the comparison with direct detection experiments we found that favorable region of 
light DM is also allowed in (^4 model. The large part of the model parameter space is very 
constrained while two other allowed areas are masses above order of 165 GeV and masses around 
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Higgs resonance. The indirect detection experiments give little less stringent bounds allowing the 
masses above around 25 GeV. 

Quaternion model can, like any scalar DM model, be further investigated via invisibly decaying 
Higgs boson in vector boson fusion (VBF) and ZH production channels at ATLAS experiment 
|33j . There is a possibility of searching for monochromatic gamma ray lines from annihilation of 
particles via virtual SM Higgs into two gamma ray lines [S]. 

We can conclude that allowed discrete non-Abelian group obtained from breaking of continuous 
SU(2) gauge symmetry can lead to new dark matter phenomena and new observations. Findings 
through that kind of new phenomena make steps toward unmasking DM nature. 
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